AD-A228  000 


DTIC 


CURVATURE  COMPUTATIONS 
ON  SURFACES  IN  n-SPACE* 

Jung-Hong  Chuangf 
and 

Christoph  M.  Hoffmann:}: 

Computer  Sciences  Department 
Purdue  University 
Technical  Report  CSD-TR-1015 
CAPO  Report  CER-90-34 
September,  1990 


pit:?  v? i  no n_st * tt  v. nrr  a 

Approved  ivj  p-.-ri’c 

Djr.  ir  UpJhr.rt<H 


*  This  research  was  supported,  through  the  Leonardo  Fibonacci  Institute,  by  the  Instituto  Tren- 
tino  di  Cultura,  Trento,  Italy. 

f  Department  of  Computer  Science  and  Information  Engineering,  National  Chiao  Tung  Univer¬ 
sity,  Hsinchu,  Taiwan,  Republic  of  China.  Supported  by  the  National  Science  Foundation  under 
grant  CCR  86-19817,  and  by  the  Office  of  Naval  Research  under  contracts  N00014-86-K-0465 
and  N00014-90-J-1599. 

$  Department  of  Computer  Science,  Purdue  University,  West  Lafayette,  IN  47907.  Supported 
by  the  National  Science  Foundation  under  grant  CCR  86-19817,  and  by  the  Office  of  Naval 
Research  under  contracts  N00014-86-K-0465  and  N00014-90-J-1599. 


90  ID  H 


Curvature  Computations  on  Surfaces  in 

n-Space* 


« 


Jung-Hong  Chuang*  Christoph  M.  Hoffmamr 

August  1990 


Abstract 

Let  /bea  two-dimensional  manifold  in  n-dimensional  space,  ai.u 
let  tt(  J’')  be  its  projection  into  the  subspace  of  three  of  the  variables  in 
which  T  has  been  expressed.  We  give  an  algorithm  that  computes  the 
normal  curvature  of  ^(T)  directly  from  the  equations  of  T  without 
variable  elimination.  We  also  comment  on  applications  in  computer- 
aided  geometric  design  (CAGD). 


,  -  1  ■ 

1  Introduction 

Many  surface  operations  in  CAGD^derive  new  surfaces  from  given  ones  sub¬ 
ject  to  certain  constraints.  Examples  include  offset  surfaces,  where  a  distance 
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constraint  must  be  observed,  and  spherical  blends,  where  a  curvature  and  a 
continuity  constraint  must  be  satisfied.  While  there  exist  intuitive  descrip¬ 
tions  of  the  resulting  surfaces  that  are  easily  grasped,  a  precise  mathematical 
representation  of  the  surfaces  appears  difficult  to  obtain  in  practice. 

Indeed,  closed-form  representations,  are  available  in  principle  with  the 
help  of  elimination  theory  —  e.g.,  [4,  5]  —  or  using  Grobner  bases  tech¬ 
niques  —  e.g.,  [2.  3.  14].  In  practice,  however,  closed  forms  are  usually 
unobtainable  because  the  elimination  problems  that  must  be  solved  are  well 
beyond  the  capabilities  of  machines  and  algorithms  available  to  date:  see. 
e.g.,  [16].  In  consequence,  surface  operations  including  offsets  and  blends 
have  been  treated  individually  in  the  literature,  and  specific  approximation^ 
methods  for  the  resulting  surfaces  have  been  given  that  are  not  general/' For 
example,  offset  computations  are  addressed  in  [7,  10.  12.  13.  18,  19,  22.  25, 
26],  and  many  ingenious  and  useful  techniques  for  analyzing  and  approximat¬ 
ing  offsets  have  been  derived.  But  the  methods  proposed  in  those  papers  for 
deriving  and  analyzing  offsets  do  not  apply  unchanged  to.  say,  the  treatment 
of  other  surfaces  whose  definition  also  involves  distance  constraints,  such  as 
Voronoi  surfaces  [6,  11]. 

In  [6,  8,  11.  15,  16]  we  have  given  a  uniform  method  for  deriving  an  exact 
representation  of  offsets,  blends,  equal-distance  surfaces,  and  so  on.  The 
representation  is  a  system  F  =  0  of  nonlinear  equations  in  n  variables,  where 
n  >  3.  with  the  property  that  the  surface  of  interest  is  the  natural  projection 
of  the  solution  set  of  the  system  into  the  subspace  spanned  by  the  first  three 
variables.  Thus,  the  surface  is  conceptualized  as  the  projection  of  a  certain 
2-surface  in  n-space.  and,  as  discussed  in  [6,  15.  16],  the  auxiliary  variables  in 
the  system  F  —  0  have  a  concrete  geometric  meaning  that  can  be  exploited 
for  instance  in  engineering  design  applications. 

In  [14,  15],  a  uniform  method  for  evaluating  the  intersection  of  such  sur¬ 
faces  has  been  described.  In  this  paper,  we  develop  a  uniform  method  for 
determining  the  surface  curvature  at  a  given  point.  That  is,  we  present  an 
algorithm  for  the  following  problem:  A  surface  T  is  given  as  a  system  of  m 


o 


nonlinear  equations  in  n  variables 

—  xn)  =  0 

f 2 \£\'  ■ 

...Xn)  =  0 

/mi*  U- 

...Xn)  =  0 

It  is  assumed  that  the  solution  set  of  the  system  is  locally  a  smooth  2-manifold 
in  Rn.  and  so  we  have  ordinarily  m  =  n  —  2;  however,  as  discussed  in  [16],  in 
certain  situations  m  >  n  —  2  is  desirable.  Let  tt(  F)  be  the  natural  projection 
of  F  into  the  (x(,  x2,  x3)-subspace.  p  a  point  on  F.  and  v  a  tangent  direction 
to  F  at  the  point  p.  Determine  the  normal  curvature  of  z(F)  at  zip)  in  the 
direction  rr(  v). 

If  the  nonlinear  equations  fk  are  algebraic,  our  problem  can  be  solved  in 

principle  by  eliminating  x4 . xn  from  the  equation  system,  followed  by  the 

well-known  curvature  computation  in  3-space  using  the  shape  operator,  see. 
e.g.,  [21].  Such  a  solution  would  not  be  practical,  however,  for  in  CAGD 
applications  the  elimination  of  the  auxiliary  variables  x4 . xn  usually  can¬ 

not  be  carried  out  in  practice.  Hence,  an  algorithm  is  needed  that  avoids 
elimination  altogether.  Such  an  algorithm  is  given  here. 


2  An  Example  from  Surface  Blending 

To  illustrate  how  surfaces  can  be  described  as  sets  of  nonlinear  equations,  we 
derive  a  ruled  surface  K  that  arises  in  certain  approaches  to  blending  two 
surfaces,  e.g.,  [20,  23].  The  approach  can  be  conceptualized  as  follows:  We 
are  given  two  surfaces  (or  patches  of  surfaces)  /  and  g,  on  which  two  link 
curves  C\  and  C2  have  been  specified.  Note  that  /  and  g  are  usually  assumed 
to  be  parametric,  but  this  assumption  is  inessential  and  can  be  dropped.  The 
blending  surface  h  should  be  tangent  to  /  in  the  curve  C\,  and  tangent  to  g 
in  the  curve  C2.  See  Figure  1  for  an  illustration. 

Now.  the  approach  is  to  put  the  points  of  C\  into  i-1  correspondence 
with  the  points  of  C2,  and  to  connect  corresponding  points  px  and  P2  with  a 
straight  line  L.  Thereupon,  a  plane  P  through  L  is  considered  that  intersects 
/  and  g  in  two  curves.  Cj  and  C0 ,  that  are  blended  as  curves  by  a  curve  C\ 


Figure  1:  Blending  Surface  and  Link  Curves 


that  depends  on  L  and  on  the  way  in  which  P  has  been  chosen.  Note  that 
Ch  can  be  obtained  as  a  functional  blend  [20 J  or  as  a  parametric  curve  [23]. 
See  also  Figure  2.  The  plane  P  can  be  suitably  chosen,  say  by  requiring  that 
it  is  normal  to  /  at  p1(  as  illustrated  in  Figure  3.  The  approach  effectively 
reduces  the  three-dimensional  surface-blending  problem  to  a  two-dimensional 
curve-blending  problem,  for  it  can  then  be  proved  that  the  surface  obtained 
in  this  manner  is  a  blending  surface  for  /  and  g.  Moreover.  [21]  proves  that 
curvature  continuity  of  the  curve  blend  C\  implies  curvature  continuity  of  the 
surface  h  with  /.  provided  that  the  tangents  to  C\  and  C\  do  not  coincide  at 


Figure  2:  Blending  Curve  C\  in  the  Plane  P 


Figure  3:  The  Plane  P  through  the  Line  L 


px.  That  is.  the  planes  P  must  not  contain  the  tangents  of  the  link  curves. 

A  difficulty  with  this  approach  to  constructing  a  blending  surface  is  to 
find  a  simple  method  for  establishing  the  correspondence  of  points  on  the 
link  curves  Ct  and  C2.  In  [23].  Pegna  proposes  the  following  idea:  Design  a 
space  curve  Co.  say  as  a  Bezier  curve,  and  let  the  link  curves  be  its  orthogonal 
projection  onto  the  two  surfaces.  That  is.  C\  is  the  orthogonal  projection  ol 
Co  onto  /.  whereas  C>  is  the  orthogonal  projection  of  Co  onto  g.  Then  two 
points  pi  and  p2  on  Ct  and  C2  correspond  precisely  when  they  are  the  image 
of  the  same  point  p0  on  Co-  See  also  Figure  4. 

The  lines  L  under  this  point  correspondence  define  a  ruled  surface  A  .  and 
we  describe  a  representation  of  I\  as  a  system  of  nonlinear  equations.  This 
is  an  example  of  the  surfaces  considered  in  this  paper,  and  the  derivation 
illustrates  our  methodology  for  deriving  complex  surface  representations  by 
devising  systems  of  equations  with  auxiliary  variables.  Note  that  the  method¬ 
ology  is  akin  to  the  equational  programming  paradigm  in  programming  lan¬ 
guages.  e.g.,  [IT],  except  that  here  the  equations  must  be  interpreted  over  a 
field  whereas  in  programming  language  research  the  equations  are  typically 
understood  over  a  free  algebra. 

We  begin  the  surface  derivation  by  assigning  variable  names  to  generic 


Figure  4:  Point  Correspondence  on  C,  and  C2 
point  coordinates.  Let 


Po  =  (u0,  Co,  U’o) 

be  a  generic  point  on  Co, 

Pi  =  (Ul.t’i.U'i) 

the  projection  of  p0  onto  /. 

P2  =  (tl2,  v  2,  IV  >  ) 

the  projection  of  p0  onto  g.  and  let 

q  =  (x.y,c) 

be  a  generic  point  on  the  line  L 

through  p  1  and  p2- 

Assume  further  that  the  curve  Co  is  given  parametrically  with  the  coordinate 
lunctions 

C0:  {  H\(r),  II2(r),  Hj(r)) 

The  surfaces  /  and  g  mav  be  given  parametrically  or  implicitly.  For  the  sake 
of  illustrating  both  cases,  we  will  assume  that  /  is  given  implicitly,  as 

f{x,y,z)  =  0 

and  that  g  is  given  parametrically,  by 

g  :  (gi{s.t).g2(s.t),g3{s.t)) 

It  will  then  be  clear  how  to  modify  the  equations  in  case  both  surfaces  are 
parametric  or  both  are  implicit. 
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We  now  compile  the  equations  defining  the  system  K  by  translating  all 
geometric  constraints  on  the  points  into  equations.  For  each  constraint,  we 
obtain  one  or  more  equations  as  follows: 


The  point  p0  is  on  Co- 


uo  =  Hx(r) 
vo  =  H7(r) 
w0  =  F/3(r) 


The  point  p1  is  on  f: 
The  point  p2  is  on  g: 


The  point  q  lies  on  the  line  L. 
and  L  contains  p,  and  p2: 


Pi  is  the  projection  of  p0 
onto  /: 


p2  is  the  projection  of  pfJ 
onto  g: 


/(“  i.L'i.  W\  )  =  0 
u2  =  Pits. t) 

'■>  =  </2Us.  n 
tm  =  J/3JS.  n 


.r  =  Au,  -j~  ( 1  —  A)«2 

y  =  At*i  +  ( l  —  A)u2 
-  =  Auq  +  ( 1  —  A)ll>2 

(u0-  u I,  Vo  -  w0-  Wi)  •  (0. |£)  =0 

(u0  -  t’0  -  rt.u?o  -  «’i)  •  _Tr)  =  0 

(«0  -  ux,  c0  -  e,.  ti>o  -  ivi)  ■  (-|£,  5^7^°^  =  0 


i  no  —  ~  r: *  /r- ) 


(  U 0  —  "  •>.  c0  -  c  >.  U’o  —  J 


dn\ 

>)  o? 

«-*T» 

)s 

-*S  1 

ds 

'72 

Oqy 

)t  ' 

■M 

So.  the  ruled  surface  K  with  points  q  =  [x,y,z\  as  been  represented  by  15 
equations  in  the  16  variables 

x,  y,  z,  uo,  Vq.  Wq.  t/j .  i'| ,  w i .  «2,  c2.  ic2.  r .  s.  f.  A 
For  an  example,  see  the  appendix. 

There  is  a  redundancy  in  the  equation  system  that  has  been  introduced  on 
purpose,  when  expressing  the  constraint  that  p{  is  the  orthogonal  projection 


onto  the  implicit  surface  f.  Here 


fci  = 

t3  =  (~|^>S7>0) 


are  three  tangent  vectors  to  /  at  the  point  pi,  which  is  evident  when  the 
inner  product  with  the  gradient  vector 

Vf  =  (  — .  — ,  —L) 

dux  Ov ,  ’  0w{ 

is  computed.  Clearly,  there  can  be  only  two  linearly  independent  vectors 
among  the  t^,  but  since  some  of  the  partial  derivatives  could  vanish  at  px,  we 
cannot  decide  a- priori  which  ones  are  independent.  Instead,  it  is  convenient 
to  have  an  algebraic  dependence  in  the  system  and  adjust  the  algorithms 
that  work  with  the  system  accordingly.  See  also  [16]. 


3  Curvature  Computations 

We  assume  given  a  point  p  =  (pi,...,pn)  that  satisfies  the  system  F  —  0  of 
m  equations  in  n  variables: 


/  i(*i . xn)  =  0 

f2{x  i,...,xn)  =  0 

fm  (Xi,...,Xn)  —  0 

Moreover,  we  assume  that  the  hypersurfaces  /*,.  =  0  are  smooth  at  p,  and 
that  the  Jacobian  of  the  system. 


££l 

i IL  .. 

. 

\ 

)X2 

9ln 

*h. 

2A  .. 

.  9h 

dx\ 

3X2 

Jin 

dfm 

3fm.  .  . 

l  frn 

/ 

dx\ 

3x2 

din 
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has  rank  n  —  2.  so  that  the  solution  set  defined  by  F  —  0  is  locally  a  two- 
dimensional  manifold. 


Let  7r  be  the  natural  projection  function  from  Rn  to  R3  that  maps  points 
q  =  (xt,....£n)  to  points  7r (q)  =  (x1,x2l-E3)-  If  F  is  the  solution  of  the 
system  F  —  0  in  Rn  and  x(F)  its  projection  into  R3.  then  we  will  show 
that  the  normal  curvature  of  n(F)  at  7r(p)  is  determined  by  a  certain  linear 
combination  of  the  second  fundamental  forms  of  the  hypersurfaces  fk,  where 
0  <  k  <  m.  Moreover,  we  will  show  how  to  compute  this  linear  combination. 
In  consequence,  we  can  determine  at  ir(p)  all  curvature  properties  of  F) 
without  the  need  to  determine  7r t  F)  explicitly.  The  following  lemma  explains 
how  to  determine  the  normal  to  n{F)  at  the  point  7r(p),  in  R3. 

Lemma  1  Let  nt  be  the  normal  vector  to  the  hypersurface  fk  at  p  t  R", 
for  1  <  k  <  m.  Let  ak  be  such  that  the  last  n  —  3  components  of 

m 

n0  =  Q^nk  ~  (a,  6,  c,  0 . 0) 

fc=i 

are  zero.  Then  r(n0)  is  normal  to  ir(F)  at  ?r(p). 

Proof  If  t  =  (ti,  t2,  t3, ...)  is  tangent  to  T  at  p,  then  7r(t)  =  (tl,t2<h) 
is  tangent  to  n(f)  at  x(p).  Since  7r(t)  •  7r(n0)  =  0.  n0  is  normal  to  x(F)  at 

7T(p).  □ 

Meusnier’s  theorem  in  differential  geometry  states  that  the  second  deriva¬ 
tive  of  any  curve,  through  a  point  p  that  lies  on  a  smooth  surface  /  in  3-space, 
has  a  projection  onto  the  surface  normal  that  does  not  depend  on  the  curve. 
The  theorem  generalizes  to  higher  dimensions,  where  it  can  be  stated  as 
follows: 


Lemma  2  (Meusnier)  Let  3{t)  be  a  curve  on  the  hypersurface  /  in 
R"  through  the  point  p  =  d(0),  and  let  v  =  <9#(0 )/dt.  Then 


Tp(v)  •  v  =  — VviV  = 


d2d(0) 

OF 


n 


where  n  is  the  normal  to  /  at  p,  and  N  is  the  (unit)  normal  vector  field  of  /. 


Note  that  Lp(v)  is  the  shape  operator.  Therefore,  we  can  define  the 
normal  curvature  of  the  hypersurface  /  at  p  in  the  tangent  direction  v  as 


«(v)  =  Lp(v)  •  v 


9 


where  v  has  unit  length.  We  now  develop  the  main  theorem. 

Consider  a  curve  3(t)  on  the  2-surface  F  in  Rn  that  contains  p,  ana 
assume  that 

m  =  p 

am  _ 

at  ~  v 

Then  the  curve  lies  on  every  hypersurface  fk,  for  1  <  k  <  m.  and  v  is  tangent 
to  F  and  each  /<... 

The  curve  3  projects  to 

*(3(t))  =  iJi(t),32{t)<Mt)) 


We  choose  numbers  a*,  such  that 

m 

n0  =  aknk  =  (a,  b.  c.  0 . 0) 

k=  1 


where  the  are  normals  to  fk  at  p,  and  a2  -f  b2  +  c2  —  1.  That  is.  the 
projected  normal  7r(no)  has  unit  length.  We  will  abbreviate  with  3  and 
with  1 3.  With  L(pfc)  the  shape  operator  of  fk  at  p,  we  obtain 


£r=1<H(L<fc)(v).v)  = 


zr=i  Qfc(/^(0)  ■  n*) 
Er«i(^(o)  ■  ak^k) 
3(0)-(T.^aknk) 

J(0)  •  (a.  6. c.  0 . 0) 

(A(0). /i2(0), /?3(0))  •  7r( n0) 


(1) 


So.  if  the  curve  3{t)  has  been  parameterized  such  that  the  vector 

*(v)  =  (ft(0),/?2(0),&(0)) 

has  unit  length,  then  expression  (1)  is  the  normal  curvature  of  ff(F)  at  x(p) 
in  the  direction  7r(v).  We  summarize  this  fact  as  follows: 

Theorem  3  Let  p  be  a  poin ,  on  the  2-surface  F  in  R",  where  F  is  the 
intersection  of  m  smooth  hypersurfaces  fk,  and  assume  that  F  is  smooth  at 
p.  Let  v  =  (ylt ....  vm)  be  a  tangent  vector  to  F  at  p,  with  v2  +  v\  +  v$  =  l. 
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and  let  there  be  numbers  a t,  such  that  n0  =  OfciH,  where  xi  n0)  has  unit 
length  and  its  last  n  —  3  components  are  zero.  If  L is  the  shape  operator 
of  the  hypersurface  fk  at  p.  then 

k  =  jr^kL^iv)  ■  v 

A:=l 

is  the  normal  curvature  of  ~{F)  at  n(p)  in  the  direction  m  v ) . 

In  consequence,  the  following  algorithm  computes  the  normal  cu^-  ature 
of  the  projected  surface 

1.  Adjust  the  length  of  the  tangent  vector  v  such  that  its  projection  tti  v) 
has  unit  length. 

2.  For  each  hvpersurface  j\  compute  the  normal  vector  ru  at  p. 

3.  Find  a  unit  normal  to  the  projected  surface  by  solving  the  linear  system 

entailed  by  ££1=1  it  =  (a,  6.  c,  0 . 0),  and  then  adjusting  the  a*  such 

that  a2  +  b2  +  c*  =  l. 

4.  For  1  <  k  <  m,  compute  L^’(v)  ■  v. 

5.  Compute  k  =  (£™=l  a tZ^fc,(v))  ■  v. 

It  is  straightforward  to  implement  this  method.  A  different  algorithm  for 
computing  the  curvature  of  the  projected  surface  can  be  given  also,  based  on 
finding  a  parametric  curve  J(t)  on  T,  see  [8]. 

Note  that  the  computation  of  L(k\\)  ■  v  can  be  based  on  the  following 
observation. 

Lemma  4  Let  g  =  0  be  a  hypersurface,  and  let  Ar  =  Vp  be  the  normal 
vector  field  of  g  at  the  point  p.  With  H  the  Hessian  matrix  of  g  at  p,  if 
v  =  (t?!, ....  vn)  is  a  tangent  vector  of  g  =  0  at  p,  then  Lp{\)  ■  v  =  -vTHv 

Proof 

£p(v)  •  v  =  —  VyiV  •  V 

=  -  Vy  Vp  ■  v 
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-  V 


On 


V\ 


On 


Ox  t  Ox  2 

On 


..  V, 


On 

0  X  ri 


[  \0x  1 

n  ;y2  _  n 

E-^»>vE 


v-  , 

=  -ZoTfir1”*'*' 

i . ;  =  L  T  ^ 

-  -vr//v 


(£) 

(p)  ■  V.. 

-r(£) 

ip) 

o-i > 

(p)c,.  .  . 

n  ■  \2 

■(  pic 

Ox,Ox  2 

’  TH  0x,0xn 

Because  of  the  bilinearity  of  the  form  —vrHv.  notice  that  wecan  rephrase 
I  heoretn  3  in  a  form  better  suited  to  computing  the  normal  curvature  at  tt( p) 
in  different  directions. 

Corollary  5  Let  p  be  a  point  on  the  2-surface  T  in  Rn.  where  J-  is  the 
intersection  of  rn  smooth  hypersurfaces  /*,  and  assume  that  JF  is  smooth  at  p. 
Let  v  =  i  r,, ....  cm)  be  a  tangent  vector  to  T  at  p,  with  rf  +  rt  +  =  1.  and 

let  there  be  numbers  or-  such  that  n0  =  satisfies  the  hypotheses 

of  Theorem  3.  With  n<,.  =  V fk  and  IIK  the  Hessian  of  the  hypersurface  /&  at 
p.  let 

//»  =  T  a,//, 

t=i 

Then 

K  =  -Vr//0V 

is  the  normal  curvature  of  ~(J-)  at  ~{p)  in  the  direction  tt\  v). 

4  Normal  Curvature  of  a  Parametric  Sur¬ 
face 

We  illustrate  our  results  by  deriving  a  formula  for  the  normal  direction  and  for 
the  normal  curvature  of  parametric  surfaces,  in  the  direction  of  the  isopara¬ 
metric  lines.  Another  example  is  given  in  the  appendix. 
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We  consider  a  parametric  surface 

x  =  h  1  ( s ,  t ) 
y  =  h2(s.t) 

=  =  h3{s,t) 

as  the  projection,  into  ix,  y.  z  )-space.  of  the  2-surface  in  5-dimensional  lx.  ?/,  c.  ,■ 
space  obtained  by  intersecting  the  three  hvpersurfaces 

/i  :  x  —  /( l(s,  f ) 

h  :  y-hlisA)  (2) 

/a  :  -  -  ht(s.t) 

Ihe  coordinate  tunctions  hi.  h2.  h'i  are  assumed  to  be  analytic  and  twice 
differentiable.  This  assumption  holds  in  particular  for  the  rational  poly¬ 
nomial  functions  used  in  CAGD.  We  abbreviate  the  partial  derivatives  of 
the  coordinate  functions  by  subscripts:  for  example,  we  write  h 2,  instead  of 
0h2/ds. 

The  normals  to  the  j\  are 


V/x  =  (1,0,0, -hl„-hlt) 

V/2  =  (0, 1,0,  —h2a,  -h2t)  (3) 

V/3  =  (0,0, 1,  —  /i  3 , ,  — /i3() 

Applying  Lemma  1.  we  must  solve  a  linear  system  in  order  to  obtain  the 
normal  direction  to  the  parametric  surface: 

ctj/il,  T  ct2h2,  -j-  0:3/185  =  0  . 

(%\h\t  4-  a2h2t  -p  a3/i3(  —  0 

The  system  is  solved  by 

at  =  h2,h3t  —  h3,h2t 

a2  =  hZ,hlt  ~  1  (5) 

o3  -  hl,h2t  —  l>2, hlt 

In  consequence,  the  normal  vector  to  the  parametric  surface  is  given  by 

(01,02,03)  =  (/tl„/i2„ /i3s)  x  (hltJi2t.  h3t) 
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To  obtain  the  unit  normal  n0  in  projection,  we  must  adjust  the  by  dividing 
by  the  length  of  the  cross  product,  i.e.. 


no  =  \  hl„  h'2,,  h3s)  x  (h\t,  !i2>.  h'it)/m 

m  =  i|(/il,,  Ii2,,  hi,)  x  1 7tlt. /t2,. /i3()|| 

Next,  we  compute  the  normal  curvature  in  the  tangent  direction  of  the 
isoparametric  line  t  =  const.  The  Hessians  of  the  system  (2)  are  the  matrices 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  -/il„  [st 

0  0  0  -hl,t  -h  1„ 

0  0  0  0  0  \ 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  — /i2„  -h:2,t 

0  0  0  -A‘2„  / 

0  0  0  0  0  \ 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  -/*3„  -hiH 

0  0  0  -h3,t  -h:]tt  ) 

Therefore,  the  matrix  H o  is  siven  by 

/  0  0  0  0  0  \ 

0  0  0  0  0 

H0  =  0  0  0  0  0  (6) 

0  0  0  a  b 
\  0  0  0  6  c  ) 

where 

a  =  -(Al,f,/i2„,/i3„)  ■  (M„/i2„/i3,)  x  (hlt,h2uhZt)/m 
b  =  -(hl,tJi2,tJi3,t)  ■  (hl„h2„hZs)  x  (hlt,h2„hZt)/m 
c  =  -{h\ttJi2n,h3tt)  ■  (h\„h2,Ji3,)  x  hZt)/m 

m  =  \\{hl,Ji2,Jil,)  x  (hlt,h2t,h3t)\\ 
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Now  the  vector 


v  =  {hl,Ji2,.h3,,  1.0) 

is  tangent  to  the  2-surface  defined  by  (2)  and  projects  to  the  tangent  of  the 
isoparametric  line  t  =  const.  When  divided  by  n  =  j(/il,.  /i2,.  /i3,)||,  the 
projected  vector  has  unit  length.  In  consequence,  the  curvature  is  given  by 

^"f=const  —  v  H0vjn 
=  — a/n 2 

Similarly,  the  normal  curvature  in  the  direction  of  the  isoparametric  line 
.-i  =  const  is 

,  > 

^  <»= const  '  ~~~f  fTl 


5  Summary 

We  have  presented  a  method  for  determining  the  local  curvature  of  the  pro¬ 
jection.  into  3-space,  of  a  2-surface  in  n-space.  The  method  does  not  rely 
on  expensive  elimination  computations,  and  we  have  implemented  it.  Our 
techniques  are  useful  in  situations  in  which  complex  surfaces  cannot  be  rep¬ 
resented  in  a  simple  closed  form,  or  are  easily  approximated  by  parametric 
surfaces,  when  the  surfaces  instead  have  been  expressed  using  systems  of 
equations  in  more  than  three  variables.  Examples  of  such  surfaces  include 
olfset  surfaces,  Voronoi  surfaces,  fixed  and  variable-radius  spherical  blending 
surfaces,  and  auxiliary  surfaces  such  as  the  ruled  surface  A  of  Section  2  or 
the  trimming  surfaces  (11]  used  in  the  definition  of  the  skeleton.  Algorithms 
such  as  the  one  presented  here  are  part  of  an  infrastructure  of  surface  in¬ 
terrogation  methods  that  should  make  surface  representations  by  systems  of 
equations  a  reasonable  alternative  allowing  practical  work  with  geometrically 
constrained  surfaces. 
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Figure  5:  The  surfaces  /  and  g.  the  curve  C0,  a  generator  of  the  ruled  surface, 
and  corresponding  points. 

Appendix:  Curvature  of  a  Ruled  Surface 

We  give  an  example  of  a  ruled  surface  of  the  type  described  in  Section  2  and 
compute  its  curvature  in  the  direction  of  the  ruling.  Note  that  the  curvature 
should  then  be  zero. 

Let  /  be  a  cylinder  of  radius  1  whose  axis  is  parallel  to  the  z-axis  through 
the  point  f  — 1.0.0).  Let  g  be  a  cylinder  of  radius  3  whose  axis  is  the  y- 
axis.  We  assume  that  /  is  given  implicitly  and  that  g  is  given  parametrically. 
Furthermore,  let  Co  be  a  circle  in  the  plane  z  =  4  of  radius  3/2  centered 
at  (  —  1/2. 0.4).  C0  is  also  given  parametrically.  We  formulate  the  equations 
of  Section  2  that  describe  the  ruled  surface  that  passes  through  the  two 
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orthosonai  projections  oi  C  o  onto  /  and  g. 


u0 

>'o 

(Co 

(  U  \  -j-  1  —  1 

111 
I'l 
(C  2 
X 

!J 

2(w0  -  'C,)(u  1  +  1) 
—  2(  ti o  —  "i  )0  +  2( vo  ~  '’i  )(ui  +  U 

2(u>0  -  (c,)i’i 

—  12s(  Ui)  —  ui)  +  6(1—  •->')(  Wo  —  n' 2 ) 

(’o  ~  r2 


3(1  -  c~  i  _  !_ 

2(l  +  r-|  2 

3r 

1  +  r' 

1 

0 

3(1  -  >-*» 

1  + 

/ 

(is 

1  +  s’ 

\ ii  i  +  (  1  -  A  iu2 

Af]  -t-  ( 1  —  A )  t’2 
Aa’i  +  ( 1  —  A  la’2 

0 

0 

0 

0 

0 


Note  that  there  are  16  variables  and  15  equations. 

On  the  surface  defined  by  the  equations,  we  choose  the  point 


p  —  ( y,  c.  u(),  t'o,  te'o,  ut,  ci,  aq,  U2,  t’2,  tt’2,  r,  s.  A)  _ 

3  ....  6  ......  3  „  12  „  V17  - 


_ ..  0.2  +  -+T-,  1,0.  4.0. 0.  4.  -7=t>  0.  -+^-0. 

2  V 1 7  Vl<  vli  Nli 


-O.i) 


The  curve  and  surfaces  are  shown  in  two  projections  in  Figure  5.  The  points 
in  the  figure  are  the  following  projections  of  p: 


q  =  {x,y,z) 

Po  =  (u0,  uo^’o) 

Pl  =  («i,  C'i,  u?, ) 

P2  =  (U2,V2,W2) 
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Note  that  q  —  \px  +  ( 1  —  \  ip2-  For  the  point  p  in  R16.  we  construct  the  sum 
of  the  Hessians.  The  resulting  matrix  H0  is  as  follows 
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q,  =  0.2776737  a2  =  0.303813S 

a,  =  0.2372107  a4  =  0.415S092 

os  =  0.8316184  a«  =  0.5553474 

o  7  =  0.4770799  o8  =  0.564125S 

o9  =  1.4232639  a0  =  1.7147960 

At  the  point  p.  the  tangent  vector  in  16-dimensional  space  is 
3  12 

v  =  (-=.0.  ~=  -  4.  0.0.0.0.0.0,0.0,0.0.0.-n 
n/T7  v  17 

and  projects  in  3-space  to  the  vector 

,  3  12 

( -7=. 0.  -=  -  4) 

Vll  >/l7 

which  is  in  the  direction  of  the  ruling.  We  compute 

v//0vr  =  0 

so  the  normal  curvature  of  the  projected  surface,  in  the  direction  of  the  ruling. 

is  zero. 
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